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摘要 和 氧 原子 钟 具有 较 高 的 短期 稳定 度 , 将 其 作为 主 钟 可 在 短期 内 产生 高 精度 的 本 地 时 间 信 号 . 但 氧 钟 存在 
频率 漂移 现象 ， 导 致 其 长 期 稳定 度 较 差 ， 从 而 影响 本 地 时 间 的 准确 性 . 毫秒 脉冲 星 自转 高 度 稳定 , 借助 于 其 长 
期 稳定 度 高 的 特性 ,可 定期 实现 对 氧 原子 钟 的 频率 驾驭 ， 并 对 实时 信号 加 以 控制 . 首先 分 析 了 国际 脉冲 星 计 时 
阵 (International Pulsar Timing Array, IPTA) 第 二 批发 布 数据 中 四 颗 毫 秒 脉冲 星 的 稳定 度 随时 间 的 变化 , 同时 
采用 哈达 玛 方差 分 析 了 中 国 科学 院 国家 授时 中 心 (National Time Service Center, NTSC) 一 台 氧 钟 的 频率 稳定 
PERE, 最 终 给 出 了 利用 脉冲 星 驾驭 氧 原子 钟 频率 的 方法 . 
关键 词 脉冲 星 , mus TP, 时 间 , 方法 : 数据 分 析 

中 图 分 类 号 : P127; ”文献 标识 码 : A 


1 gl 言 国际 原子 时 (International Atomic Time, TAT) 
时 间 频 率 技术 的 发 展 在 国民 经 济 、 国 防 建设 系统 是 由 全 球 共同 维护 的 时 间 基 准 系 统 , 它 是 国际 


和 基础 科学 中 发 挥 着 重要 作用 , 通信 、 互 联网 以 及 权 度 局 (Bureau International des Poids et Mesures, 
卫星 导航 等 技术 均 需 要 建立 精密 、 统 一 的 时 间 系 BIPM) 利 用 全 球 90 多 个 时 频 实验 室 的 500 多 台 原 子 


统 . 为 提供 时 间 尺 度 服务 , 美国 海军 天 文 台 (United 。” 钟 建立 和 维持 的 融 精 度 时 间 系统 . TAI 由 ALGOS 算 
States Naval Observatory, USNO), 德国 物理 技 TUIS, 输入 数据 为 全 球 实验 室 参 与 归 算 的 原子 
术 研 究 所 (Physikalisch Technische Bundesanstalt， 钟 与 UTC (PTB) 的 比 对 钟 差 (此 处 UTC (Coordi- 
PTB), 日 本 信息 与 通信 技术 研究 所 (National Insti- nated Universal Time) 为 协调 世界 时 , UTC (kR 
tute of Information and Communications Technol- 表 实 验 室 k 保 持 的 地 方 协 调 世 界 时 )， 比 对 方式 为 
ogy, NICT), 中 国 科学 院 国 家 授时 中 心 (National GPS (Global Positioning System) 共 视 技 术 、 卫 星 
Time Service Center, NTSC) 等 均 建立 了 自己 的 守 双向 时 间 频 率 传递 技术 以 及 GPS 载波 相位 时 间 传 
时 系统 . 守 时 系统 建立 的 目的 在 于 实现 统一 的 时 间 递 技术 . 加 权 平 均 算法 的 关键 是 确定 原子 钟 的 权重 ， 
频率 基准 , 并 保证 基准 的 稳定 性 和 准确 性 . 稳定 性 ALGOS 算 法 属于 典型 的 加 权 平 均 算法 , 该 算法 中 
征 的 是 振荡 器 在 一 定时 间 内 产生 同样 时 间 和 频 原子 钟 的 权重 采用 月 速率 方差 的 倒数 计算 . 同时 ， 
率 的 能 力 , 而 准确 性 表征 的 是 与 理想 值 的 关系 , 时 为 限制 单 台 钟 权 值 过 大 对 原子 时 标 造 成 不 利 的 影 
间 基 准 的 建立 需 同时 兼顾 准确 度 与 稳定 度 吓 . 响 , ALGOS 算 法 采用 了 最 大 权 约 束 原则 . 经 加 权 平 
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均 计 算 后 得 到 的 原子 时 间 尺 度 噪 声 水 平 会 有 效 降 
低 ， Bese BOER 高 于 单 台 原子 钟 . 

实验 室 的 原子 钟 包 括 守 时 型 原子 钟 和 基准 型 
原子 钟 . 目前 参与 UTC 计算 的 守 时 钟 主 要 为 氨 钟 和 
4& p. 为 更 好 地 发 挥 氧 原子 钟 优 良 的 短期 稳定 度 和 
低 噪声 特性 ， 从 2014 年 开始 ，BIPM 降 低 了 钨 钟 参 
与 TAI 计 算 的 总 权重 , BI TAPES, 截至 
2021 年 12 月 , 氧 钟 在 计算 自由 时 标 (Echelle Atom- 
ique Libre, EAL) 中 所 占 的 权重 达到 了 85% 左 右 , 成 
为 UTC 保 持 的 核心 守 时 钟 . 由 于 氢 钟 的 精度 不 断 提 
高 , 绝 大 多 数 实验 室 的 主 钟 频率 源 由 钨 钟 更 换 为 氢 
钟 . 2005 年 3 月 起 , 国家 授时 中 心 采用 氧 钟 做 主 钟 . 
对 氧 钟 开 展 相应 的 研究 , 有 望 进一步 提高 本 地 时 间 
尺度 的 性 能 
地 球 时 (Terrestrial Time, TT) 作 为 一 种 理想 的 


= 


时 间 尺 度 , 其 时 间 单 位 与 大 地 水 准 面 上 的 国际 单位 
M (Système International d'Unités, SI) £5 — 8, Jr 


点 定义 为 在 1977 年 1 月 1 日 零 时 (TAD, TT = TAI + 
32.184 s. TAI 为 TT 的 一 种 实现 ,考虑 到 实时 性 需 
求 ，TAI 每 月 归 算 一 次 ， 因 此 TT (TAIT) 的 精度 在 
一 定 程度 上 受到 影响 . 且 由 于 主 频 标的 不 确定 度 
以 及 频率 鸭 驭 过 程 的 不 完美 , TAI 存在 长 期 的 频率 
波动 . 此 外 ，TAI 还 存在 季节 性 波动 向. 为 克服 TT 
(TAD 的 上 述 缺 陷 , BIPM 利 用 了 所 有 可 用 的 基准 钟 
的 数据 , 在 计算 过 程 中 进行 了 黑体 辐射 频 移 改 正 ， 
并 对 EAL 与 TT 间 的 频率 偏差 进行 低 通 滤波 与 插值 
处 理 , 扣除 了 EAL 频 率 的 周年 项 与 长 期 变化 项 , 最 
终 得 到 TT 的 最 优 实现 一 TT (BIPM)! TT (BIP- 
MXX) 每 年 归 算 一 次 , XX 为 年 份 的 后 两 位 数字 . 
TT (BIPM) 的 钟 差 数据 发 布 周期 为 一 年 ， 影 
响 时 间 比 对 的 实时 性 , 而 脉冲 星 有 望 在 缺失 比 对 
数据 的 情况 下 保证 本 地 时 间 的 准确 性 . 脉冲 星 自 
转 非常 稳定 ， 其 旋转 产生 的 周期 性 脉冲 信号 可 被 
地 面 射电 望远镜 接收 . 这 种 稳定 的 周期 信号 可 作 
为 一 种 天 然 的 频率 源 ， 称 为 脉冲 星 钟 56. 毫秒 脉 
冲 星 自转 的 长 期 稳定 度 可 媲美 高 精度 的 原子 钟 ， 
利用 参考 钟 对 所 有 脉冲 星 计 时 残 差 影响 相同 的 特 
性 , 可 构建 综合 脉冲 星 时 间 尺 度 (Ensemble Pulsar 


Timescale, EPT)'® 9，EPT 是 同时 观测 多 颗 毫 秒 
脉冲 星 ， 并 设计 相应 的 算法 提取 计时 残 差 的 公共 
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信号 得 到 . EPT 可 用 于 检验 原子 时 的 长 期 波动 n9， 
并 作为 TT 的 一 种 新 的 实现 方式 ， 记 为 TT (PT), 
PT (Pulsar Time) 为 脉冲 星 时 . 原子 时 的 建立 与 
水 冲 星 时 的 建立 属于 两 个 完全 不 同 的 物理 过 程 ， 
TTT (PT) 与 TT (BIPM) 的 差异 很 小 L171， 因此 
其 他 时 间 尺 度 相 对 于 TT (BIPM) 的 波动 可 利用 TT 
(PT) 加 以 检验 . 
本 地 时 间 UTC (k) 是 基于 单 台 物 理 钟 产生 的 时 
间 信 号 , 而 任何 一 台 物 理 钟 相对 于 标准 时 间 尺 度 均 
td ie Es] IE ig Eh FE Bh HE TT LETS C 目前 ， 
各 实验 室 利 用 BIPM 前 期 发 布 的 UTC (k ) 与 UTC 的 
时 刻 差 , 推算 出 相位 微调 和 频率 微调 的 数据 , 进而 
对 UTC (k) 加 Bet 但 BIPM 发 布 数据 的 滞后 将 
对 UTC (k) 的 准确 性 产生 影响 . 作为 TT 的 另 一 种 实 
BU, 我 们 尝试 用 TT ree 59 驭 主 钟 频率 . 
为 了 监测 和 改善 伽利略 系统 时 间 的 长 期 稳定 
性 ,2018 年 12 月 , 欧 空 局 位 于 荷兰 的 技术 中 心 开 始 
运行 基于 脉冲 星 时 斩 驭 伽利略 卫星 导航 系统 主 控 
站 主 钟 的 “Pulchron”(Pulsar Chron)Ji H U91, 基于 
此 项 目 与 五 百 米 口径 球面 射电 望远镜 FAST (Five- 
hundred-meter Aperture Spherical radio Tele- 
scope) 等 国内 望远镜 的 高 灵敏 度 观 测 ， 周 庆 勇 等 
人 给 出 了 可 提高 北斗 卫星 导航 系统 时 间 基 准 长 期 
稳定 性 的 中 国 脉冲 星 时 地 面 服务 系统 (China Pul- 
sar Time Ground Service System, CPTGSS) 的 构 
EDRN., Hp corp E KE FY E ENTE 
阵列 开展 , 国际 上 主要 的 脉冲 星 计时 阵 包 括 澳 大 利 
亚 Parkes 脉 冲 星 计 时 阵 PPTA (Parkes Pulsar Tim- 
北美 纳 赫 兹 引力 波 计时 阵 NANOGrav 
(North American Nanohertz Observatory for 
Gravitational Waves) 以 及 欧洲 脉冲 星 计时 阵 EP- 
TA (European Pulsar Timing Array), 三 者 共同 
组 成 了 国际 脉冲 星 计时 阵 IPTA (International Pul- 
sar Timing Array). IPTA 的 主要 科学 目标 为 探测 
其 低频 引力 波 以 及 建立 脉冲 星 时 间 尺 度 ， 其 对 毫 
秒 脉冲 星 脉冲 到 达 时 间 TOA (Time of Arrival) 的 
典型 测量 精度 达到 几 百 纳 秒 . 本 文 从 IPTA 第 二 批 
发 布 数 据 IPTA2 中 选取 T ASTE ESS i BE BO ea H3] 2e 
PD fk HE (J0437-4715, J1713--0747. J1744-1134. 
J1909-3744) 216.2 yr 的 实测 数据 , 对 四 颗 星 的 稳定 
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度 进行 分 析 . 选取 的 研究 对 象 为 NT 


韩 孟 纳 等 : 脉冲 


SC 的 一 台 主 动 


型 氧 脉 泽 钟 (H1404926), 钟 差 数 据 来 源 于 BIPM 网 
站 1. 基于 此 钟 的 钟 差 数 据 , 本 文 给 出 了 利用 毫秒 脉 
冲 星 观测 数据 每 三 个 月 、 每 半年 和 一 年 为 周期 轰 
驭 一 次 主 钟 频率 的 方法 , 并 对 四 颗 星 的 综合 驾驭 能 
力 进 行 评 佑 . 


2 ”四 颗 脉 冲 星 稳定 度 评估 


本 文 


首先 对 IPTA2 的 实测 数据 [9 进行 处 理 , 得 


涉及 一 系列 的 


到 四 颗 星 的 拟 合 前 计时 残 差 , 该 过 程 


时 间 尺 度 转换 . 实际 记录 TOA 的 钟 


Tj. 为 扣除 参考 钟 对 脉冲 星 计时 精度 的 影 


钟 的 读数 需 转 到 TT (BIPM) 时 间 尺 


扣除 几何 时 延 、 光 线 弯曲 、 大 和 气 传播 时 延 、 


效应 以 及 狭义 相对 论 的 “ 动 钟 变 慢 ” 


为 测 站 的 本 
fl, ASH 
EA. 此 外 , 为 
视差 
、 广 义 相 对 论 


te E 


FET UEP Ph BEES A BUSCAS 


中 的 “引力 红 移 ”对 不 同位 置 钟 的 走 ; 
TOA 


centric 
型 


心 SSB 


f&, 


行 模 型 
差 如 图 
Zub 


致 ， 对 观测 数据 i 


时 残 差 为 有 


zn 


p" 


最 终 转 到 太阳 系 质心 坐标 时 TCB 


影响 等 因 


A, 
(Bary- 


Coordinate Time)!6 171, 脉冲 星 的 计时 模 


(Solar System Barycenter) 


参数 拟 合 , 最 终 得 到 四 颗 星 
1 所 示 , 图 中 MJD (Modified 
KH. 为 与 下 文 氧 钟 钟 差 数据 的 时 间 
行 截取 , 截取 后 数据 的 时 


!4 参 数 通过 拟 合计 时 残 差 获 得 , 可 用 于 在 太阳 系 质 
处 预报 TOA. 计 
永 冲 实际 到 达 时 间 与 模型 预报 时 间 的 差 
般 在 TCB 时 间 尺 度 下 计生 
冲 星 数据 处 理 软件 0 对 四 颗 星 


. 采用 Tempo2 脉 
的 观测 数据 重新 进 
GWEF Wh 
Julian Date) 7j 


区 间 一 


间 跨 度 


为 MJD: 53739 一 56014. 从 计时 残 差 结果 来 看 ， 除 
脉冲 星 J0437-4715, 其 他 三 颗 星 的 计时 残 差 均 在 堆 
附近 波动 , 说 明 这 三 颗 星 不 含 明显 的 计时 噪声 . 


4r 1 60r 
—— J0437-4715 
3r 4 40 L 
2r I | 20r 
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2 aA Ly La. | 
L i | We M1 bi rers B 
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oO 2f 
— 1 1 1 1 1 -60 1 1 1 1 1 
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= 30 T MJD T T T MJD, T 
O) —— J1744-1134 10 - 
£E 2. ] 
E 5. 
vor] 10r 
+ 
uc ' THN | ] 0r 
zo] sell 
O aol 中 5r 
[au . 
-20 L -10 + 
一 一 J1909-3744 
-30 + -15 上 
53500 54000 54500 55000 55500 56000 53500 54000 54500 55000 55500 56000 
MJD MJD 
于 1 四 颗 脉 冲 星 的 拟 合 后 计时 残 差 图 
Fig.1 The post-fit timing residuals of four pulsars 


*https://webtai.bipm.org/ftp/pub/tai/data/. 
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脉冲 星 的 稳定 度 采用 Matsakis 等 人 提出 的 统 定 度 结果 进行 分 析 , 可 以 看 出 脉冲 星 的 稳定 度 随 着 
计量 os 表征 L191, 其 不 要 求 数 据 的 等 间隔 分 布 , 适用 时 间 增 加 , 即 其 具有 优良 的 长 期 稳定 度 . 四 颗 星 整 
于 脉冲 星 观测 数据 的 分 析 . cs: 按照 下 式 计算 : 体 稳定 度 均 大 臻 服从 lgo。 x 5lgr (u = 一 3), 上 
E » 述 关系 对 应 图 2 中 斜率 为 -1.5 的 直线 , 说 明 四 颗 脉 

nt ie (1) ” 冲 星 的 主要 噪声 类 型 为 白 噪 声 . 对 于 观测 时 间 小 于 

一 年 的 数据 , PSR J1744-4715 的 稳定 度 略 逊色 于 其 

其 中 ,cs 为 对 长 度 为 + 的 观测 数据 进行 三 次 多 项 式 他 三 颗 星 . DOK MP Be (30437-4715. J17134-0747. 


函数 拟 合 得 到 的 三 次 项 系数 ; 
权 平 均 , 权重 为 拟 合 值 cs 不 确 
上 式 分 别 计算 脉 ; 


式 中 的 尖 括 号 表示 加 


eer 


J1744-1134. J1909-3744)256 yr 观测 数据 的 稳 


c 


AC IX 


上 定 度 的 平方 倒数 . OK 分 别 为 7.95 x 10-19. 4.38 x 10716, 4.94 x 10-15、 
星 的 稳定 度 , 得 到 四 颗 星 的 稳 3.66 x 10-19. 随 着 观测 设备 灵敏 度 的 增加 , 四 颗 星 


定 度 随 时 间 的 变化 关系 如 图 2 所 示 . 对 四 颗 星 的 稳 观测 数据 的 稳定 度 均 得 到 进一步 的 提高 . 
49710 ; 10710 r 
- =- 0437-4715 aa - =~-J17134+0747 
Ail 10 E 
10 
-12 | 
一 101? F 一 
B È gal 
N N 10 
> 49°13; 5 
i04 
10714: 10775 "Se 
10715 1 i i 
107 10"! 10° 10! 107 10"! 10° 10! 
7lyt 7lyt 
ig? 
1071! L 
-12. 
zo 
N 
S 1078 | 
104: 
10715 1 1 i i 
10? 107 109 10! 10? 10"! 10° 10! 
yr Tjyr 
图 2 ”四 颗 脉 冲 星 的 稳定 度 随时 间 的 变化 图 . 图 中 颜色 不 同 的 直线 斜率 不 同 , 蓝 色 、 红 色 、 青 色 、 洋 红色 线 的 斜率 分 别 为 2、 一 1.5、 一 1、 一 0.5. 


Fig.2 Variation of the stability of four pulsars with time. The slopes of the blue, red, blue-green and purple lines are —2, —1.5, 


—1 and —0.5, respectively. 


3 脉冲 星 驾驭 原子 钟 


NTSCRJSF BT 4 
4 FY) Ru UL RS E BE LF $68. 
Tidi AEH A LR 4 


动 型 
准确 


有 望 提 


高 时 间 


尺度 的 长 期 稳定 度 . 时 间 实 验 室 负责 


组 包括 主动 型 氧 钟 和 钨 钟 , E 产生 和 保持 的 系统 时 间 包 括 物 理 时 间 和 纸 面 时 间 . 
若 能 利用 脉冲 星 对 于 物理 时 间 , 可 利用 毫秒 脉冲 星 的 计时 信息 控制 
高 差 及 频率 漂移 并 将 其 扣除 ， 主 钟 频率 , 实时 输出 物理 信号 ; 纸 面 时 间 需 要 利用 
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历史 数据 重新 进行 原子 时 的 计算 , 以 达到 更 高 的 稳 
EE, 属于 事后 计算 . 本 文 主要 针对 主 钟 的 实时 物 
时 信号 驾 驱 开展 相应 的 研究 . 

为 提高 实时 UTC (k) 信 号 的 准确 性 , 可 利用 守 
时 钟 组 、 快 速 UTC (UTCr) 与 UTC 分 别 实现 对 主 
钟 频率 的 每 日 、 每 周 及 每 月 驾驭 R021. 为 提高 所 
钟 长 期 保持 的 准确 性 ,本文 提 出 在 上 述 驾驭 方法 
的 基础 上 , 增加 利用 毫秒 脉冲 星 的 观测 数据 定期 芍 
驭 氧 原子 钟 的 方案 . 下 文 首先 研究 了 NTSC 一 台 参 
与 TAI 归 算 的 氢 钟 H1404926 (HH 代表 原子 钟 类 型 为 
氧 钟 ，1404926 为 钟 编 号 ) 的 稳定 度 , 然后 给 出 利用 
脉冲 星 驾 驭 氧 钟 的 方法 . 


3.1 和 氢 钟 稳定 度 分 析 

氧 钟 具 有 优良 的 短期 稳定 度 , 充分 发 挥 氧 钟 的 
特性 对 于 地 方 协 SEN dA " :的 高 精度 保持 具 
重要 意义 . 图 3 给 出 了 利用 哈达 玛 方差 计算 的 
A P'H1404926 8] fa xE B£, 从 图 中 可 以 看 出 氢 钟 的 短 
稳 高 而 长 稳 较 差 , 这 主要 由 氧 钟 的 老化 引起 . 哈达 
玛 方差 采用 了 频率 数据 的 二 阶 差 分 、 相 位 数据 的 
一 阶 差 分 , 其 对 调频 闪烁 游 走 噪声 和 调频 随机 奔走 
噪声 收敛 , 且 对 于 线性 频率 漂移 不 灵敏 , 因此 适用 
于 氢 钟 的 稳定 度 分 析 . 针对 相位 数据 哈达 玛 方差 的 
定义 为 : 


a = 
N-3 : 
GXKN —3) 3 (Ging — 925,2 + 3Ti41 — vi), 
(2) 
其 中 , mu. mia. miae. Tirs AVE AMA, TH 


邻 数 据 间 的 时 间 间 隔 为 MN 为 相位 数据 的 取样 数 . 
通常 使 用 哈达 玛 方差 的 平方 根 , 即 哈达 玛 偏差 表征 
氧 钟 的 稳定 度 . 

记 单 台 原 子 钟 实现 的 TT 为 TT (Y), RAYA 
原子 钟 的 编号 . 将 本 文 使 用 的 氧 钟 实现 的 TT 记 为 
TT (H1404926), BI TT (H1404926) = H1404926 + 
32.184 s. 利用 BIPM 发 布 数据 , 经 下 式 转 换 可 以 得 
到 真实 钟 差 TT (H1404926) — TT (BIPM), 


[TT(H1404926) — TT(BIPM)] = 


— ([TT(BIPM) 
[UTC(k) 

钟 差 结果 如 图 4 所 示 ， 从 图 中 可 以 看 出 二 者 之 
间 存 在 频率 漂移 ， 可 见 单 台 钟 实现 的 TT 是 不 准确 
的 . 


10712[ 一 一 一 T — -T rr T ——] 


— TT(TAD] + [UTC — UTC(K)] + 
— H1404926]) . (3) 


H^ 
1014 rani " i at " " EM 
10! 10? 10° 
Ald 


图 3 ”哈达 玛 偏差 表征 的 氢 钟 稳定 度 随 时 间 的 变化 医 


Fig.3 The variation of stability of the hydrogen master with 


time characterized by Hadamard deviation 
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44 TT (H1404926) 5'TT (BIPM) 之 间 的 真实 钟 差 


Fig.4 The real clock difference between TT (H1404926) and 
TT (BIPM) 


Hifi PRA RI, 氧 钟 存在 频率 漂移 , 氧 钟 的 预报 
模型 采用 二 次 三 项 式 P3, 即 : 


h(t) = A+ B(t — to) 


P stri 
os 
c 
© 
wa 
N 
Pene. 
A 
x 
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KE, A(t) At Zl zt Pp EST Zl FR A, ANto IN ZI 其 中 Lp = 1.550519768 x 10-8. MERTA, & 
钟 差 ， 巨 、C 分 别 为 钟 的 频率 偏差 以 及 频率 漂移 ， ”将 TCB 时 间 尺 度 下 氢 钟 的 频率 偏差 记 为 /， 即 TCB 
对 TT (BIPM) 与 氧 钟 之 间 的 真实 钟 差 每 月 进行 二 (H1404926) — TCB (PT) 的 斜率 , 则 在 TT 时 间 尺 
次 三 项 式 拟 合 , 得 到 氢 钟 的 频率 偏差 不 是 常数 , 而 “， 度 下 氧 钟 的 频率 偏差 为 !(1 — Lp). 该 速率 差 远 小 于 
是 随时 间 变 化 的 . 脉冲 星 对 氧 钟 频率 偏差 的 估计 误差 10-15， 因 此 对 
国家 授时 中 心 采用 类 ALGOS 算 法 计算 本 地 时 ”合计 结果 的 影响 可 以 忽略 . 
间 ， 多 台 原 子 钟 的 读数 经 加 权 平 均 算 法 得 到 自 1 TOALITT (H1404926) 为 参考 时 , 得 到 四 颗 星 
时 标 后 , 结合 BIPM 网 站 发 布 的 钟 差 数据 计算 驾驭 ”的 拟 合 前 计时 残 差 如 图 5 所 示 ， 脉冲 星 数据 长 度 约 


E 


E, 对 主 钟 频 率 进行 轨 驭 , 脉冲 星 对 氢 钟 的 鸭 驭 只 62. 其 变化 趋势 基本 与 图 4 中 真实 钟 差 的 趋势 一 
需 在 记录 TOA 的 参考 钟 与 守 时 原子 钟 之 间 建 立时 致 , 证 实 了 脉冲 星 驾 驭 原子 钟 的 可 行 性 . 从 图 5 可 以 
间 比 对 网 , 也 可 通过 参考 共同 的 时 标 , 如 北斗 导航 ” 看 出 , 以 同一 台 氧 钟 为 参考 钟 记录 TOA 时 , 四 颗 肪 
卫星 的 系统 时 间 (BeiDou Time, BDT), 4 —42*.— 冲 星 的 拟 合 前 计时 残 差 变化 趋势 一 致 ; 但 由 于 每 颗 
立 联系 . 此 外 , 可 将 脉冲 是 实现 的 TT (PT) 作 为 TT — 星 本 身 的 计时 噪声 不 同 , 因此 各 星 计时 残 差 中 随机 
(BIPM) 的 一 种 检验 方式 , 并 在 BIPM 的 比 对 数据 更 项 的 波动 程度 不 同 ， 下 文采 用 抗 差 估计 [627 降低 
新 不 及 时 的 情况 下 维持 本 地 时 间 的 准确 性 . 下 面 我 。 ”脉冲 星 计 时 噪声 对 驾驶 结果 的 影响 . 

y yyh E yh ee Arh ZS Ey] SZAN NA FE 

ANA FD ct ESTE Vb} = 4S. RAE EA 324 单 星 抗 差 估计 结果 


周期 进行 驾驭. 

将 脉冲 星 时 作为 参考 时 标 , 通过 估计 原子 钟 相 
3.2 脉冲 星 钟 对 氢 钟 的 频率 轰 驭 方法 对 于 脉冲 星 时 的 频率 , 对 原子 钟 的 频率 进行 驾驭 ， 
通过 设计 参考 钟 与 主 钟 系统 频率 源 之 间 的 比 。 最终 将 得 到 一 个 稳定 性 更 高 的 时 间 尺 度 . 对 于 以 氧 
对 链 路 , 可 实现 脉冲 星 对 主 钟 频率 的 驾驶 , 以 此 提 ” ” 钟 作为 记录 脉冲 到 达 时 刻 的 参考 钟 时 得 到 的 拟 合 
高 本 地 时 间 系 统 的 准确 性 . 若 参 考 钟 无 钟 差 、 模 型 AIRE, 利用 抗 差 估计 [39 每 三 个 月 为 一 个 计 
参数 准确 且 脉 冲 星 不 存在 计时 噪声 , 则 拟 合 前 计时 算 周 期 按 (4) 式 和 欠 代 估计 参考 钟 的 频 偏 与 频 漂 , 直至 
残 差 将 呈现 出 均值 为 零 的 白 噪声 特征 . 为 研究 脉冲 结果 趋 于 稳定 . 在 采用 加 权 平 均 算法 计算 原子 时 的 
星 对 氧 钟 的 驾驶 结果 , 本 文 将 TOA 的 参考 时 间 尺 度 过 程 中 , 原子 钟 的 频 偏 与 频 漂 估 计 一 般 采 用 最 小 二 
转换 为 TT (H1404926). 由 于 计时 残 差 是 在 TCB 乒 “。 乘 算 法 . 但 最 小 二 乘 算法 受 粗 差 影响 较 大 , 若 测量 
度 下 计算 得 到 的 , 因此 ， 拟 合 前 计时 残 差 应 为 TT 数据 中 混入 少量 粗 差 , 则 会 导致 估计 结果 与 真 值 的 
(H1404926) 相 对 于 TT (PT) 的 变化 在 TCB 时 间 尺 偏离 程度 较 大 . 抗 差 估计 是 在 观测 值 超过 正常 范围 


度 下 的 对 应 值 .为 得 到 TT (H1404926) ~ TT (PT), 。 时 ， 限 制 观测 值 对 平 差 的 影响 . 采用 抗 差 估计 可 以 
需 计算 TT 与 TOB 之 闻 的 比例 关系 , TOSH OE 在 一 定 程度 上 减 小 脉冲 星 本 身 的 计时 噪声 对 估计 


标 时 TCG (Geocentric Coordinate Time) 2198095 。 结果 的 影响 , 提高 频率 偏差 与 频率 漂移 的 估计 精度 


异 包括 周期 项 与 长 期 漂移 项 , 即 ”: 抗 差 估计 的 权 因 子 函数 应 用 较 广泛 的 有 IGG3 函 数 
TCB- TCG =- to (TT, +TT,)+A,, (5) Huber 89), IGG3K BN: 
LB 
其 中 Lc = 1480826867 x 1078, TT 为 线性 项 , TT, 1 TE 
Y 3 , JI X ? 

为 非 线性 项 ，A; 为 观测 者 依赖 项 . TT 与 TCG 的 束 NET 
率 之 比 为 一 常数 ， ppt z1- Lol. TT5ETCB does t TRU ) y^ C1 < |u| Sc, (7) 
的 走时 速率 存在 如 下 关系 : os 

dTT 

arog ~!~ 4s: (6) (7 =1,2,--+,n). 
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45 ”参考 钟 为 TT (H1404926) 时 四 颗 星 的 拟 合 前 计时 残 差 氏 
Fig.5 The pre-fit timing residuals of four pulsars when the reference clock is TT (H1404926) 
Huber Zt 7j: 图 6 与 图 7 所 示 , 其 中 图 6 使 用 的 权 因 子 函 数 为 IJGG3 
函数 , 图 7 使 用 的 权 因子 函数 为 Huber 函 数 . 从 图 6 与 
ej < e; 图 7 的 对 比 结果 来 看 , 二 者 对 应 的 频率 偏差 与 频率 
w= i Tl |v;| >c, (8) 漂移 的 趋势 基本 一 致 , 说 明 两 种 权 因 子 函 数 都 可 以 
a 有 效 抑制 脉冲 星 噪 声 对 估计 结果 的 影响 . 为 方便 起 
bd Jh, 下 文 我 们 只 讨论 使 用 IGG3 函 数 的 结果 . 从 四 颗 
AF, ;为 每 个 TOA 的 序号 , wj 为 权 因 子 , vj 为 观测 星 的 估计 结果 来 看 , J1744-1134 与 其 他 三 颗 星 的 估 
值 与 拟 合 值 之 差 , 即 残 差 . 若 以 正 态 分 布 统计 量 构 计 结 果 偏 离 程度 较 大 . 这 是 因为 该 星 的 计时 精度 相 
造 权 函数 ， 则 在 特定 的 显著 性 水 平 下 权 函 数 的 界 ”” 较 于 其 他 三 颗 星 较 低 . 因此 , 在 脉冲 星 驾 驭 原子 钟 
值 可 以 国定 . 为 尽 可 能 地 减 小 脉冲 星 噪声 对 结果 的 的 应 用 过 程 中 , 应 选择 计时 精度 高 且 不 受 红 噪 声 影 
影响 并 保留 原子 钟 本 身 的 特性 , 根据 拟 合 后 残 差 wj 响 的 毫秒 脉冲 星 . 
的 区 间 范 围 选取 合适 的 界 值 . 对 于 J1744-1134 脉 冲 3.2.2 ”四 颗 星 加 权 平 均 结果 
星 , IGG3 函 数 的 误差 控制 界 值 cu 、cz 分 别 选 为 10-8 为 了 有 效 降 低 单 颗 星 的 噪声 对 估计 结果 的 影 
51075, Huber Zt f] Z8 c7310-5; 对 于 其 他 三 颗 响 ， 下 面 对 四 颗 星 的 拟 合 结果 进行 加 权 平 均 处 理 ， 
E, IGG3 函 数 的 界 值 cu、c? 分 别 设置 为 5 x 10-85 每 颗 星 的 权重 反比 于 各 星 以 TT (BIPM) 为 参考 时 
1079, Huber 函 数 的 界 值 c 设 为 5 x 1078. 得 到 的 拟 合 后 计时 残 差 均 方 根 值 RMS (Root Mean 
抗 差 估计 得 到 的 频率 偏差 与 频率 漂移 结果 如 Square) 的 平方 . 进行 总 权重 归 一 化 处 理 后 , 得 到 四 


64 4 


颗 星 (J0437-4715、 — J17134-0747. J1744-1134, 
J1909-3744) HA HW; 4) ll 0.228. 0.290. 0.025. 
0.457, 结果 如 表 1 所 示 . 四 颗 星 综合 处 理 的 结果 如 
图 8 所 示 . 同时 , 图 8 还 给 出 了 每 季度 利用 UTC 估 计 
的 氢 钟 H1404926 频 偏 与 频 漂 值 以 及 TT (BIPM) 对 
应 的 估计 结果 .通过 与 TT (BIPM) TT 
(H1404926) 得 到 的 拟 合 结果 对 比 可 知 ，UTC 对 氧 
钟 的 频 偏 与 频 漂 估计 结果 与 TT (BIPM) 估 计 的 结 
果 基 本 一 致 . 而 经 四 颗 脉冲 星 观测 数据 估计 的 氨 钟 
频 偏 与 氨 钟 相对 于 TT (BIPM) 的 频 偏 在 前 期 差别 
较 大 , 后 期 符合 程度 较 好 , 这 主要 与 脉冲 星 观测 设 
备 灵敏 度 的 提升 有 关 . 脉冲 星 时 估计 的 氧 钟 频率 漂 
移 在 TT (BIPMJ) 对 应 估计 值 附 近 波 动 ， 波动 程度 
前 期 较 大 而 后 期 较 小 , 波动 量 与 脉冲 星 自 身 噪 声 、 
观测 设备 噪声 、 传 播 路 径 导 致 的 噪声 以 及 脉冲 星 
计时 模型 精度 有 关 . 为 了 与 每 季度 驾驭 结果 进行 
比较 , 图 9 与 图 10 分 别 给 出 了 每 半年 与 一 年 利用 TT 
(BIPM)、PT 及 UTC 对 氧 钟 H1404926 的 频率 偏差 、 
频率 漂移 的 估计 结果 . 从 估计 结果 来 看 , 驾驭 周期 
越 长 , 脉冲 星 时 的 驾驶 结果 与 TT (BIPM) 的 驾驭 结 
果 越 接近 . 为 了 更 好 地 比较 估计 结果 , BOT 
周期 为 一 年 时 三 者 的 估计 结果 展示 在 表 2 中 , KP 
第 一 列 为 每 个 驾驭 周期 的 起 始 时 间 , B. CARI 
表 频 率 偏 差 与 频率 漂移 的 估计 值 . 从 表 2 中 的 结果 
可 以 看 出 , 三 者 的 频率 偏差 估计 结果 相近 . 在 后 期 
研究 中 , 我 们 将 筛选 计时 噪声 小 的 多 颗 毫 秒 脉冲 上 
进行 综合 处 理 , 并 采用 其 他 方法 更 好 地 扣除 脉冲 星 
自身 因素 对 估计 结果 的 影响 . 


表 1 四 颗 脉 冲 星 的 归 一 化 权重 
Table 1 The normalized weights of the four 


pulsars 
PSR RMS/us Wj 
J0437-A715 0.273 0.228 
J1713+0747 0.242 0.290 
J1744-1134 0.832 0.025 
J1909-3744 0.193 0.457 


X 学 dk 6 期 
IRH EOS est Eh ATE, 可 编制 相应 


的 软件 , 将 利用 脉冲 星 时 估计 的 频率 偏差 与 频率 漂 
移 量 作为 频率 补偿 值 加 载 到 相位 微调 仪 上 , 实现 对 
主 钟 频 率 的 定期 驾驭 , 从 而 提高 本 地 时 间 的 准确 度 . 
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(b) The estimated frequency drift 
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Fig.6 The frequency deviation and frequency drift of the 
hydrogen maser H1404926 obtained by robust estimation 
using data of the four pulsars. The weight factor function is 
IGG3 function. 
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Fig.7 The frequency deviation and frequency drift of the hydrogen maser obtained by robust estimation using data of the four 


pulsars. The weight factor function is taken to be the Huber function. 
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图 8 四 颗 脉 冲 星 估计 结果 综合 处 理 得 到 的 氧 钟 频率 偏差 与 频率 漂移 . 图 中 的 黑色 数据 点 为 由 TT (BIPM)-H1404926 每 三 个 月 进行 二 次 三 项 式 拟 
合 得 到 的 拟 合 值 , 蓝 色 数据 点 为 UTC-H1404926 每 三 个 月 的 估计 结果 , 红色 数据 点 为 四 颗 脉 冲 星 加 权 平 均 得 到 的 估计 值 . 


Fig.8 The frequency deviation and frequency drift of the hydrogen maser obtained by comprehensively processing the 
estimation results of the four pulsars. The black data points in the figure are the values obtained by the quadratic trinomial 
fitting of TT (BIPM)-H1404926 every three months, the blue data points are the estimated results of UTC-H1404926 every 


three months, and the red data points are the estimates obtained by the weighted average of the results of the four pulsars. 
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Fig.9 The frequency deviation and drift of the hydrogen maser obtained by comprehensively processing the estimation results 


of the four pulsars. The steering cycle is half a year. 
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Fig.10 The frequenc 
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y deviation and drift of the hydrogen maser obtained by comprehensively processing the estimation results 


of the four pulsars. The steering cycle is one year. 
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2 驾驭 周期 为 一 年 时 频率 偏差 与 频率 漂移 的 估计 结 


Table 2 The estimation of frequency deviation and drift when sone cycle is one year 


B 
MJD 


C/Hz 


T (BIPM) UTC PT 


T (BIPM) UTC PT 


53739 4.9186 x 10? 4.9498 x 10713 5.2578 x 10713 
54104 —2.6852 x 10 !? —2.4909 x 1071? —2.5066 x 10718 
54469 —1.9812 x 1071? —1.9827 x 1071? —2.0649 x 107"? 
54834 —2.9659 x 10 !? —2.9470 x 1071? —2.9978 x 10^? 
55199 —3.8617 x 1071? —3.8572 x 10? —3.8804 x 10? 
55564 —4.7142 x 10? —4.6960 x 1071? —4.6948 x 10^? 


4 总 结 与 展望 

本 文 利 用 脉冲 星 长 期 稳定 度 高 的 特性 , 给 出 了 
脉冲 星 时 驾驭 单 台 氧 钟 的 方法 . 利用 IPTA 第 二 批 
发 布 数据 中 的 四 颗 脉 冲 星 数 据 , 根据 文中 方法 得 到 
氨 钟 的 频率 偏差 与 频率 漂移 . 此 外 , 采用 加 权 平 均 
算法 分 别 得 到 了 以 每 季度 、 半 年 、 一 年 为 周期 脉 
冲 星 对 氢 钟 的 综合 驾驭 结果 , 并 与 TT (BIPM)RS AS 
驭 结果 进行 对 比 . 由 于 脉冲 星 后 期 数据 TOA 测 量 精 
度 高 于 前 期 , 二 者 的 频率 偏差 估计 结果 在 所 用 数据 
段 的 后 期 基本 一 致 , 说 明 脉冲 星 源 的 选择 、 脉 冲 星 
的 观测 精度 以 及 数据 处 理 方法 等 均 会 对 驾驭 结果 
产生 影响 . 此 外 ,对 不 同 驾驭 周期 的 结果 进行 对 比 
发 现 ， 轰 驭 周期 越 长 PT 与 TT (BIPM) 的 驾驭 结果 
接近 程度 越 高 , 说 明 脉冲 星 的 优势 主要 体现 在 其 长 
期 稳定 度 . 本 文 结果 证 实 脉冲 星 驾 驭 所 钟 具 有 可 行 
TE. 该 方法 可 以 在 一 定 程度 上 弥补 TT (BIPM) 发 布 
周期 长 的 不 足 ， 有 望 提高 本 地 时 间 UTC (k) 的 准确 
TE, 并 对 后 期 的 相关 研究 具有 一 定 的 参考 价值 . 
脉冲 星 TOA 的 测量 精度 及 脉冲 星 钟 模型 参数 
情 度 是 影响 脉冲 星 驾 驭 原子 钟 精度 的 重要 因素 . 为 
提高 脉冲 星 时 轰 驭 原子 钟 的 水 平 ， 需 将 计时 精度 较 
高 的 毫秒 脉冲 星 作为 候选 源 , 同时 提高 观测 设备 的 
性 能 ， 以 提高 TOA 的 测量 精度 . 此 外 ,， 需 对 候选 脉 
冲 星 开 展 密集 而 持续 的 观测 , 以 缩短 TOA 的 采样 间 
隔 , 并 实现 在 时 间 上 对 原子 钟 的 持续 驾驭 . 后 续 工 
ER, 脉冲 星 筛选 、 候 选 星 取 权 方法 以 及 将 脉冲 星 
时 与 原子 时 融合 构建 综合 时 间 尺 度 的 算法 也 是 值 


ns 
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—2.3454 x 1077 —2.3462 x 10-?9? —2.4878 x 10-79 


—4.0756 x 1077 —4.0715 x 10-20 —4.2095 x 10-?9 


—2.9962 x 10-?? —2.9679 x 107° —2.7114 x 107?° 


—3.0414 x 10-?? —3.0438 x 10-?9? —2.9494 x 10-?9 


—2.6645 x 1077 —2.6570 x 107° —2.5691 x 10-?? 


—2.6781 x 10-?? —2.6838 x 10-?9? —2.7819 x 10-79 


得 研究 的 课题 . 
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Steering the Frequency of Hydrogen Master Based on Pulsar 
Observations 
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(2 Key Laboratory of Time and Frequency Primary Standards, Chinese Academy of Sciences, Xi’an 710600) 
(3 School of Astronomy and Space Science, University of Chinese Academy of Sciences, Beijing 100049) 


AnsBsrRACT Hydrogen masters have high short-term stability, using them as master clocks can generate 
high-precision local standard time in a short period of time. However, the hydrogen maser has a frequency 
drift phenomenon, which leads to its poor long-term stability, thus affecting the accuracy of the local time. 
The rotation of the millisecond pulsar is highly stable. With the help of its high long-term stability, it 
can regularly control the frequency of the hydrogen atomic clock, thus control the real-time signal. The 
stability of four millisecond pulsars in the second data set released by the International Pulsar Timing 
Array (IPTA) is analyzed. At the same time, the frequency stability of a hydrogen maser of the National 
Time Service Center (NTSC) is analyzed by Hadamard variance. Finally, a strategy for steering the 
frequency of hydrogen master using pulsars is given. 


Key words pulsars, hydrogen master, time, methods: data analysis 
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